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Abstract 

We  address  optimal  placement  of  vehicles  with  simple  motion,  to  intercept  a  mobile  target  that  arrives  stochastically  on 
a  line  segment.  The  optimality  of  vehicle  placement  is  measured  through  a  cost  function  associated  with  intercepting  the 
target.  With  a  single  vehicle,  we  assume  that  the  target  either  moves  with  fixed  speed  and  in  a  fixed  direction  or  moves 
to  maximize  the  vertical  height  or  intercept  time.  We  show  that  each  of  the  corresponding  cost  functions  is  convex,  has 
smooth  gradient  and  has  a  unique  minimizing  location,  and  so  the  optimal  vehicle  placement  is  obtained  by  any  standard 
gradient-based  optimization  technique.  With  multiple  vehicles,  we  assume  that  the  target  moves  with  fixed  speed  and 
in  fixed  direction.  We  present  a  discrete  time  partitioning  and  gradient-based  algorithm,  and  characterize  conditions 
under  which  the  algorithm  asymptotically  leads  the  vehicles  to  a  set  of  critical  configurations  of  the  cost  function. 

Key  words:  Coverage  control.  Pursuit-evasion  games,  Optimization,  Autonomous  and  Multi-agent  systems. 


1.  Introduction 

Vehicle  placement  to  provide  optimal  coverage  has  re¬ 
ceived  lot  of  recent  attention  due  to  potential  applications 
in  environmental  monitoring,  patrolling  a  piece  of  land  be¬ 
tween  nations  and  even  in  popular  games  such  as  soccer. 
This  work  addresses  vehicle  placement  to  minimize  a  cost 
associated  with  intercepting  a  mobile  target  that  appears 
randomly  on  a  segment. 

In  static  environments,  vehicle  placement  problems  are 
analogous  to  geometric  location  problems,  wherein  given  a 
set  of  static  points,  the  goal  is  to  find  supply  locations  that 
minimize  a  cost  function  of  the  distance  from  each  point 
to  its  nearest  supply  location  (cf.  Megiddo  and  Supowit 
(1984)  and  Zemel  (1984)).  For  a  single  vehicle,  the  av¬ 
erage  distance  to  a  random  point,  generated  according  to 
a  probability  density  function  is  given  by  the  Weber  or 
the  continuous  1  median  function,  for  which  there  exists  a 
global  minimizer  as  shown  in  Fekete  et  al.  (2005),  termed 
as  the  median.  For  multiple  distinct  vehicle  locations,  the 
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expected  distance  between  a  random  point  generated  ac¬ 
cording  to  a  probability  density  and  one  of  the  locations 
is  known  in  literature  as  the  continuous  Weber  or  the  con¬ 
tinuous  multi-median  function,  e.g.,  see  Drezner  (1995). 
For  more  than  one  location,  the  multi-median  function  is 
non-convex,  and  thus  determining  locations  that  minimize 
the  multi-median  function  is  hard  in  the  general  case.  It 
is  of  interest  to  characterize  the  set  of  critical  points  of  the 
multi-median  function.  Cortes  et  al.  (2004)  have  character¬ 
ized  the  set  of  critical  points  for  the  problem  of  deploying 
a  group  of  robots  in  a  region  to  optimize  a  multi-median 
cost  function.  Schwager  et  al.  (2009)  provide  an  adaptive 
control  law  to  enable  robots  to  approximate  the  density 
function  from  sensor  measurements.  Martinez  and  Bullo 
(2006)  presented  motion  coordination  algorithms  to  steer 
a  mobile  sensor  network  to  an  optimal  placement.  More 
recently,  Kwok  and  Martinez  (2010)  presented  a  coverage 
algorithm  for  vehicles  in  a  river  environment. 

In  mobile  target  scenarios,  the  cost  for  the  vehicle  is 
a  function  of  relative  locations,  speeds  and  motion  con¬ 
straints  considered.  For  an  adversarial  target,  the  optimal 
vehicle  motion  is  obtained  by  solving  a  min-max  pursuit- 
evasion  game,  in  which  the  target  seeks  to  maximize  while 
the  vehicle  seeks  to  minimize  a  certain  cost  function.  The 
vehicle  strategy  is  a  version  of  the  classic  proportional  nav¬ 
igation  guidance  law  (cf.  Guelman  (1971)).  With  con¬ 
straints  such  as  a  wall  in  the  playing  space  or  non-zero 
capture  distance,  strategies  with  optimal  intercept  time 
have  been  derived  in  Isaacs  (1965)  and  in  Pachter  (1987). 

We  consider  a  line  segment  on  which  a  mobile  target 
appears  via  a  known  spatial  probability  density  and  one 
or  multiple  vehicles  seek  to  intercept  it.  The  goal  is  to  de¬ 
termine  vehicle  placements  that  minimize  a  cost  function 
associated  with  the  target  motion.  This  work  is  an  exten- 
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sion  of  Bopardikar  et  al.  (2010),  where  we  introduced  the 
placement  problem  for  target  motion  with  fixed  speed  and 
in  fixed  direction,  and  for  a  uniform  target  arrival  density. 
We  address  single  and  multiple  vehicle  scenarios.  With  a 
single  vehicle,  we  consider  a  class  of  cost  functions  and  es¬ 
tablish  properties  such  as  convexity,  smoothness  and  the 
existence  and  uniqueness  of  a  globally  minimizing  vehi¬ 
cle  location.  We  show  that  the  cost  functions  associated 
with  the  target  moving  with  fixed  speed  and  in  a  fixed 
direction,  and  with  the  target  seeking  to  maximize  the 
distance  from  the  segment,  fall  in  the  class  of  cost  func¬ 
tions  that  we  have  analyzed.  The  cost  function  for  target 
motion  that  maximizes  the  intercept  time  is  shown  to  be 
proportional  to  the  continuous  1-median  function.  With 
multiple  vehicles  and  the  target  moving  with  fixed  speed 
and  in  a  fixed  direction,  we  first  provide  an  algorithm  to 
partition  the  line  segment  among  the  vehicles  and  charac¬ 
terize  its  properties.  With  the  expected  intercept  time  as 
the  cost  function,  we  propose  a  Lloyd  descent  algorithm 
in  which  every  vehicle  computes  its  partition  and  moves 
along  the  gradient  of  the  expected  time  computed  over  its 
partition.  We  characterize  conditions  under  which  the  ve¬ 
hicles  asymptotically  reach  a  set  of  critical  configurations. 

This  paper  is  organized  as  follows.  The  problem  is 
formulated  in  Section  2.  Single  vehicle  scenarios  are  ad¬ 
dressed  in  Section  3.  The  multiple  vehicle  scenario  is  an¬ 
alyzed  in  Section  4. 

2.  Problem  Statement 

We  consider  vehicles  with  simple  motion  and  speed  up¬ 
per  bounded  by  unity.  A  target  arrives  at  a  random  po¬ 
sition  (a;,0)  on  the  segment  G  :=  [0,  W]  x  {0},  termed 
the  generator ,  via  a  specified  probability  density  function 
(f)  :  [0,  W]  — ►  M>o-  We  assume  that  the  density  func¬ 
tion  <f>  is  bounded,  i.e.,  there  exists  an  M  >  0  such  that 
<f>{x)  <  M,  V x  £  [0,  W],  and  cj>  is  positive  on  a  set  of  posi¬ 
tive  measure.  The  target  moves  with  bounded  speed  less 
than  that  of  the  vehicles,  and  is  intercepted  or  captured 
if  a  vehicle  and  the  target  are  at  the  same  point.  The 
goal  is  to  determine  vehicle  placements  that  minimize  a 
certain  cost  function  based  on  the  maneuvering  abilities  of 
the  target.  Specifically,  we  consider  the  following  cases. 

2.1.  Single  Vehicle  Case 

We  determine  a  location  peRx  R>o  that  minimizes 
Cexp  :  R2  ->  M  given  by 

rW 

Cexp(p)  :=  /  C(p,x)(j)(x)dx,  (1) 

Jo 

where  C  :  M2  — >  K>o  is  an  appropriately  defined  cost  of  the 
vehicle  position  p.  In  what  follows,  we  seek  to  minimize 
the  following  different  cost  functions. 

(i)  Expected  constrained  travel  time:  We  assume  that 
the  target  arriving  at  (a;,  0)  translates  in  the  positive  In¬ 
direction  with  speed  v  <  1.  From  Bopardikar  et  al.  (2010), 


the  cost  function  C  for  this  formulation  is 


T{p.x) 


v/(l  -  v2){X  ~  x)2  +  Y2 
1  —  v2 


vY 

1^2’ 


(2) 


which  is  the  time  taken  for  the  vehicle  to  intercept  the 
constrained  target. 

(ii)  Expected  vertical  height:  The  cost  function  C  for 
this  formulation  is  the  vertical  height  H(p,x)  which  the 
target  seeks  to  maximize  before  being  intercepted. 

(iii)  Expected  intercept  time:  The  cost  function  C  for 
this  formulation  is  the  intercept  time  Ti(p,  x)  which  the 
target  seeks  to  maximize. 

Explicit  formulae  for  the  quantities  H  and  Ti  are  de¬ 
rived  in  Section  3.2,  and  are  illustrated  in  Figure  1. 
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Figure  1:  Intercepting  a  target  that  seeks  to  maximize  either  the 
vertical  height  H  or  the  time  Ti  until  intercept. 


2.2.  Multiple  Vehicles  Case 

We  assume  that  the  target  translates  in  the  positive 
Y -direction  with  speed  v  <  1.  As  shown  in  Figure  2,  given 
m  >  2  vehicles  having  complete  communication,  the  goal  is 
to  determine  vehicle  locations  p j  G  [0,  W]  x  R>0,  for  every 
i  G  {l,...,m},  that  minimize  the  expected  constrained 
travel  time  given  by 


^exp(Pl)  ■  •  ■  ;  Pm) 


j-W 

/  min  T(pi,x)<j)(x)dx, 

J  0  *£{!,. ..,m} 


(3) 


where  T(pi,x)  is  given  by  Eq.  (2). 


Pi 


Figure  2:  Intercepting  a  target  having  constrained  motion. 


3.  Single  Vehicle  Scenarios 

We  first  analyze  a  class  of  cost  functions.  This  form  will 
appear  in  two  distinct  scenarios,  the  expected  constrained 
travel  time  and  the  expected  vertical  height. 


2 


3.1.  A  class  of  cost  functions 

We  assume  that  the  cost  function  is  given  by  Eq.  (1), 
where  the  function  C  has  the  form 


C{X,  E,  x )  :=  a^b(X  -  x )2  +  Y2  -  cY,  (4) 


and  a,  b ,  and  c  are  positive  constants,  with  a  >  c. 

The  partial  derivatives  of  Cexp(X,  Y)  with  respect  to 
X  and  Y  are  given  by 


dC, 


rW 


oir=°b 


(. X  —  x)<f>(x) 

<0  a/ b{X  -  x)2  +  Y2 


dx, 


(5) 


da 


rW 


BY 


=  aY 


<t>{x) 


'o  y/ b(X-x)2  +  Y 2 


dx  —  c. 


(6) 


Lemma  3.3  (Uniqueness)  There  exists  a  unique  vehi¬ 
cle  location  (X*,Y*)  £  ] 0,  lb[xl>o  that  minimizes  Cexp. 

Proof:  Let  there  be  two  locations  (Xi,  Ei)  and  (X2,  Y2) 
that  minimize  the  expected  cost.  From  Lemma  3.1,  since 
the  expected  cost  Cexp  is  convex  in  X  and  Y,  a  con¬ 
vex  combination  of  (Xi,Ei)  and  ( X2,Y2 )  also  minimizes 
the  expected  time.  Thus,  the  necessary  conditions  for 
minimum  are  satisfied  by  (X(a),Y(a))  :=  (aXi  +  (1  — 
a)X2,  aY\  +  (1  —  o)E2),  for  every  a  £  [0, 1].  Thus, 


{X(a)-x)<j){x)  ^  =  n 

\J  (bX(a)  —  x)2  +  E  (a)2 

_ Y{a)(j)(x) _ ^  =  c 

\Jb{X(a)  —  x)2  +  Y (a)2  *  «' 


Lemma  3.1  (Convexity  of  expected  cost)  The  expected  Since  the  above  conditions  hold  for  every  a  £  [0,1],  the 


cost  Cexp(X,  Y)  is  convex  in  X  and  Y  over  the  domain 
M  x  M>o  - 


partial  derivatives  of  the  above  conditions  evaluated  at  a  = 
0,  must  equal  zero.  Thus,  upon  simplifying,  we  obtain, 


The  proof  involves  showing  that  the  Hessian  matrix  of 
C  with  respect  to  X  and  Y  is  positive  semi-definite.  The 
complete  proof  is  presented  in  ?. 

Lemma  3.2  (Existence  of  Minima)  There  exists  a  ve¬ 
hicle  location  (X*,Y*)  £  ] 0,  W[xM>o  that  minimizes  Cexp. 

Proof:  CeXp  has  continuous  partial  derivatives  in  the 
domain  ]0,  W[xR>0  (cf.  Eq.s  (5),  (6)).  We  show  that  a 
minimizer  cannot  lie  on  the  boundary  of  the  region  [0,  W]  x 
R>o-  We  begin  by  showing  that  Y*  exists  and  is  finite. 
Taking  the  limit  of  Cexp(X,Y)  as  Y  — >  +oo, 


rw 


rw 


(X2  -  x)Y2(Y1  -  Y2)  -  Y2(X1  -  X2) 
Jo  ( b{X2  -  x)i  +  E,2)3/2 

{X2  -  x)Y2{X r  -  X2)  -  (W  -  Y2)(X 2  -  x)s 
(b(X2  -  x)2  +  E2)3/2 


(j>{x)dx  =  0, 


4>{x)dx  =  0, 


where  <j>{x)/(b{X2  —  x)2  +  E22)3/2  =:  f(X2,Y2,x)  is  strictly 
positive  for  Y2  >  0.  Multiplying  the  first  equation  by  (X\  — 
X2),  the  second  by  (E  —  Y2),  and  adding  the  equations, 


i-W 

/  f(X2,Y2,x)(Y2{X1-X2)-(X2-x)(Y1-Y2))2dx  =  0. 

Jo 


p  vv 

liminf  Ceyi^{X^  Y)  >  liminf(a  —  c)Y  /  cj)(x)dx  =  +00, 
Y— >-+00  Y  —>+00  J  q 

since  by  assumption,  a  >  c.  Thus,  Y*  exists  and  is  finite. 

Finally,  to  show  that  a  minimizer  lies  in  ] 0,  W[xR>o, 
we  need  to  prove  two  statements:  (a)  Y*  ^  0,  and  (b) 
X*  e]0,  W[.  To  show  (a),  Eq.  (6)  along  with  the  assump¬ 
tion  cj>( x)  <  M,  for  every  x  £  [0,  W],  yields 


d  c. 


,-w 


<MYa 


dx 


'o  \/b(X  -  x)2  +  Y2 


—  c 


< 


MYa 
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(log(W  +  \/W2  +  Y2/b)  —  log  (Y/Vb))  —  c. 


Thus,  limsupF^0+  dCexp/dY  <  —  c.  Thus,  for  Y  near  zero 
the  gradient  of  Cexp  points  in  the  negative  E-direction, 
implying  that  E*  ^  0. 

To  show  (b),  we  first  observe  that  for  a  given  E,  in  the 
limit  as  X  — >  ±oo,  Cexp  — >  +oo,  and  therefore  X*  must 
be  bounded.  Finally,  the  claim  follows  since  the  partial 
derivative  of  Cexp  with  respect  to  X  is  strictly  negative  at 
X  =  0  and  is  strictly  positive  at  X  =  W . 

These  statements  coupled  with  the  convexity  of  Cexp 
with  respect  to  X  and  E  imply  the  existence  of  a  minimizer 
in  the  region  ]0,  W[xR>0.  ■ 


Since  f(X2,  Y2,x)  >  0,  we  must  have  Y2{X\  —  X2)  —  (X2  — 
x)(Yi  —  Y2)  =  0,  for  every  x  at  which  f(X2,Y2,x)  >  0, 
which  is  feasible  only  if  Xi  —  X2  =  0  and  Y\  —  Y2  —  0.  ■ 

We  now  present  the  main  result  for  this  section. 

Theorem  3.4  (Minimizing  expected  cost)  From  an  ini¬ 
tial  location  infix  K>o  and  by  using  a  gradient  optimiza¬ 
tion  technique,  the  vehicle  reaches  the  unique  point  that 
minimizes  the  expected  cost  Cexp. 

Proof:  The  gradient  of  Cexp  with  respect  to  X  and  E 
is  a  continuous  function  of  X  and  E  in  the  region  M  x  R>o- 
The  function  Cexp  is  convex  in  X  and  E  (cf.  Lemma  3.1) 
and  has  a  unique  minimizer  in  ]  0,  Pf[xR>o  (cf.  Lemmas  3.2 
and  3.3).  Thus,  a  gradient  optimization  technique  (cf.  Boyd 
and  Vandenberghe  (2004))  leads  the  vehicle  to  the  unique 
global  minimizer  of  CEp.  ■ 

Let  V  C  [0,  W],  and  choose  <j>{x)  such  that  <j>{x)  =  0,  Vx  £ 

[0,  W\  \  V.  Then,  Theorem  3.4  yields  the  following  result. 

Corollary  3.5  (Any  subset  of  generator)  LetV  C  [0,  W] 
have  positive  measure.  Then,  from  an  initial  location  in 
[0,  W]  x  R>o,  and  following  gradient  descent  with  V  as  the 
region  of  integration,  the  vehicle  remains  inside  [0,  W]  x 
R>0  at  all  subsequent  times. 
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Theorem  3.4  answers  the  problem  of  minimizing  the 
expected  value  of  T,  given  by  Eq.  (2),  with  a  :=  1/(1  —  v2) , 
b  :=  (1— v2)  and  c  :=  v/(l  —  v2),  and  a  >  c.  In  general,  it  is 
difficult  to  provide  closed  form  expressions  for  the  vehicle 
location  that  minimizes  the  expected  time.  A  special  case 
is  described  in  Remark  3.6. 

Remark  3.6  (Equal  speeds)  In  this  case,  the  optimal 
placement  in  the  A  variable  is  at  the  centroid  of  the  dis¬ 
tribution  </>,  with  the  optimal  Y  given  by 


Figure  3:  Move  towards  top- most  strategy  for  the  evader. 


pW  I  rW 

X*  =  J  <f(x)xdx-  Y*  =  ]l  J  <Kx)(X*  -  x)2dx.  □ 

3.2.  Optimal  Placement  for  Adversarial  Target 

We  consider  two  types  of  cost  functions  that  the  evader 
tries  to  maximize;  vertical  height  and  intercept  time. 

3.2.1.  Minimizing  the  Expected  Vertical  Height 

We  first  present  the  solution  to  the  differential  game 
with  payoff  equal  to  the  vertical  height.  If  the  evader  is 
slower  than  the  pursuer,  then  the  Appolonius  circle  is  the 
boundary  of  the  set  of  all  points  which  the  evader  can 
reach  without  being  captured.  The  following  property  is 
stated  in  Isaacs  (1965). 


strategy  and  the  resulting  optimal  vertical  height  of  the  in¬ 
tercept  point  is 

H(X,  Y,  x)  = 

1  —  v2  1  —  vz 

Comparing  the  expression  for  H  given  by  Lemma  3.8  with 
the  definition  of  C  in  Eq.  (4),  we  have  a  :=  v/(l  —  v2), 
b  :=  1  and  c  :=  u2/(l  —  v2),  and  a  >  c  since  v  <  1.  Thus, 
by  applying  Theorem  3.4,  we  obtain  the  following  result. 

Theorem  3.9  (Minimizing  expected  height)  From  an 
initial  location  in  R  X  R>o,  by  using  a  gradient  optimiza¬ 
tion  technique ,  the  vehicle  reaches  the  unique  point  that 
minimizes  the  expected  height  Hexp . 


Proposition  3.7  (Appolonius  circle  during  pursuit) 

If  the  pursuer  and  the  evader  both  travel  straight  toward 
a  point  U  on  the  Appolonius  circle ,  then  any  new  such 
circle,  obtained  from  a  pair  of  simultaneous  intermediate 
positions  of  the  pursuer  and  the  evader,  is  tangent  to  the 
original  circle  at  U,  and  is  contained  in  the  original  circle. 


3.2.2.  Minimizing  the  Expected  Intercept  Time 

In  this  formulation,  we  assume  that  the  evader  is  con¬ 
strained  to  remain  above  or  on  the  A-axis.  Thus,  the 
underlying  differential  game  in  this  set  up  is  the  classic 
wall  pursuit  game,  proposed  and  solved  in  Isaacs  (1965). 
We  present  the  main  result  for  completeness. 


The  optimal  pursuit  strategy  (cf.  Isaacs  (1965))  is  to 
choose  its  velocity  vector  such  that  the  line  joining  the 
pursuer  and  the  evader  remains  parallel  at  all  times,  while 
reducing  the  distance.  So  for  optimal  placement,  it  suffices 
to  determine  the  optimal  evader  strategy.  Algorithm  1 
summarizes  our  evader  strategy,  shown  in  Figure  3. 


Lemma  3.10  (Wall  Pursuit  game)  The  evader  strat¬ 
egy  that  maximizes  the  intercept  time  is  to  move  towards 
the  furthest  point  of  the  Appolonius  circle  on  the  X-axis. 

This  optimal  evader  strategy  is  illustrated  in  Figure  4. 


Algorithm  1:  Move  towards  top-most 
Assumes:  Pursuer  at  (X,  Y).  Evader  at  (x,  0). 
i:  Compute  center  and  radius  of  the  Appolonius  circle: 


O  :=  ( Ox,Ov )  =  ( 

R  := 


x  —  v2X  —v2Y 


1  —  v* 


1  —  v2  7  1  —  v2 
V/A  -x)2  +  Y2. 


)■ 


2:  Move  towards  the  point  (Ox,Oy  +  R )  with  speed  v. 


The  following  result  is  immediate  from  Proposition  3.7. 


Now,  given  a  convex  region  Qcl  and  a  density  func¬ 
tion  if  :  Q  — >•  K>o,  the  median  (cf.  Fekete  et  al.  (2005)) 
Pmed  is  the  unique  global  minimizer  of 


Lemma  3.8  (Move  towards  top-most  is  optimal)  The 

strategy  move  towards  top-most  is  the  evader’s  optimal 


|  p  —  z\ip(z)dz. 
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We  now  present  the  main  result  of  this  section. 


Theorem  3.11  (Optimal  location  is  the  Median)  The 

median  point  of  the  region  [0,  W]  x  {0}  with  the  density 
function  (f>  uniquely  minimizes  the  expected  intercept  time. 

Proof:  Using  Lemma  3.10  and  the  Pythagoras  theo¬ 
rem,  the  intercept  time  Ti  is  given  by 


Ti  (X,Y,x) 


x  —  vX 
1  —  v 


Figure  5:  To  determine  pairwise  dominance  regions. 


where  R  is  the  radius  of  the  Appolonius  circle  drawn  at 
the  initial  instant.  Since  pursuer  placement  on  the  X-axis 
results  into  decreasing  the  intercept  time  Ti,  we  have 

,,2  i  q,,  /-w 

Tiexp(X)  =  - - 2“  /  \X  ~  x\(j>(x)dx, 

Jo 

which  is  minimized  uniquely  by  the  median  of  the  region 
[0,  W]  x  {0}  with  the  density  function  <f>.  ■ 

4.  The  Case  of  Multiple  Vehicles 

We  now  address  the  multi-vehicle  placement  problem. 

4-1.  Dominance  Region  Partition 

We  introduce  a  generator  partitioning  procedure  by 
defining  dominance  regions  between  each  pair  of  vehicles. 

Definition  4.1  (Pairwise  dominance  region)  Fori,j  £ 
{1, . . . ,  m};  the  pairwise  dominance  region  Uij  C  [0,  W]  of 
Pi  with  respect  to  p j  is  the  set  of  target  locations  for  which 
vehicle  p,  takes  lesser  time  to  intercept  the  target  than  p j : 

Uij  :=  {x  £  [0,  W]  |  T(pj,  x)  <  T(pj,x)}. 

In  what  follows,  we  describe  a  procedure  to  determine 
Ujj,  which  is  summarized  in  Algorithm  2.  Without  loss  of 
generality,  assume  that  Xj  <  Xj.  If  Yj  =  Yj,  i.e.,  the  vehi¬ 
cles  are  at  the  same  distance  from  the  generator,  then  Uij 
is  the  piece  of  G  that  lies  in  the  half-plane  that  is  formed 
by  the  perpendicular  bisector  of  the  segment  joining  pi 
and  p j  and  which  contains  Pi.  Now  if  Y)  <  Yj,  then  we 
look  for  points  (x,0)  in  G  for  which  T(pj,x)  <  T(pj,x). 
By  setting  (1  —  v 2)  =:  b,  Eq.  (2)  gives 


Algorithm  2:  Pairwise  Dominance  Region 
Assumes:  Distinct  p j  =  ( Xi,Yi ),  Pj  =  ( Xj,Yj ). 
i:  if  Yt  =  Yj,  then 

f  [0,  (Xj  +  Xj)/2],  iiXi<Xj 


2: 

3:  else 


U^  := 


[{Xi  +  Xj)/2,W\,  ifXi  >Xj 


8  :=  arctan2 (Yj  —  Yi,  Xj  —  Xj)  +  7r/2 

£12  ■=  two  roots  of  0  =  4(sin2(0)  —  v2)£2 
+4 (Yi  +  Yj)  sin (8)£  +  (Yj  +  Yj)2  -  v2\\pi  -  p: 

2/1,2  :=  (Yi  +  Yj)/ 2  +  sin(0)4>2 
if  i/i  >  0  and  y2  >  0  then 

xq2  :=  (Xi  +  Xj)/ 2  +  cos(0)4,2 


Uij  := 


[o,  w]  n  [x1,x2\, 
[0,  W]  \  ]x\,  x2[, 


if  Yi  <  Yj 

if  Y)  >  Yj. 


10 


12 


else 


k  :=  index  in  {1,  2}  for  which  >  0 
a:  :=  (Xj+Xj)/2  +  4cos (9) 


13: 


Uij  := 


[0,  W]  n  ]— oo,  x],  if  Xj  <  Xj 
[o,  w]  n  [x,  +oo  [,  ifXj>Xj 


where  £  and  8  :=  arctan2((Y,  —  Yj),  (Xj  —  Xj))  +  7t/2  are 
as  shown  in  Figure  5.  Let  £\  and  i2  be  the  roots  of  the 
above  quadratic.  Then  the  Y-coordinates  of  the  candidate 
boundary  points  A\  and  A2  are  given  by 

[yi,  2/2]T  =  [1,  1]T(>)  +Yj)/2  +  [4,  £2}T  sind. 


\J b(Xi  —  x)2  +Y2  —  vYi  <  yj b(Xj-xY+Y2-vY3 .  (7) 

On  simplifying,  one  can  show  that  Eq.  (7)  is  quadratic 
in  x  having  real  roots,  which  provides  at  most  two  points 
for  the  boundary  between  Uij  and  Ujj .  To  determine  the 
boundary  points,  consider  the  perpendicular  bisector  of 
the  segment  joining  pj  and  p j,  as  shown  in  Figure  5.  We 
look  for  points  A\  and  A2  on  this  bisector  such  that  the 
distances  of  Ai  and  A2  from  the  real  line  is  v  times  their 
respective  distances  from  the  vehicles.  This  gives  rise  to 
the  following  quadratic  equation  in  the  variable  £ 

4(sin2  8-v2)£2+4(Yi+Yj)  sin8£  =  -(Vi-|-Y)-)2+u2||pj-pJ||2, 


Now,  Ai  and  A2  are  both  boundary  points  if  and  only  if 
both  have  positive  Y-coordinates.  It  can  be  shown  that 
there  exists  at  least  one  among  them  which  has  positive 
Y-coordinate.  There  arise  two  cases: 

(i)  If  there  are  two  candidate  points  A\  and  A2  (as 
in  Figure  5),  then  we  look  at  their  corresponding  X  co¬ 
ordinates,  (xi,x2)  given  by  Step  7.  For  (x, 0)  £  Gfl 
[xi,x2]  x  {0},  we  have  T(pj,x)  <  T(pj,x),  and  thus  Uij 
is  G  fl  [xi,  x2]  x  {0}. 

(ii)  If  there  is  only  one  candidate  point  A±,  then  we  look 
at  its  X  coordinates,  Xi  given  by  Step  9.  By  assumption 
Xj  <  Xj,  and  hence  for  (x,  0)  £  G  fl  [— oo,Xi]  x  {0},  we 
have  T(pi,  x)  <  T(pj,x)  and  thus  Uij  is  Gfl[— oo,  xi]  x{0}. 

Thus,  we  have  established  the  following  property. 
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Algorithm  3:  Dominance  region 

Assumes:  Distinct  locations  {pi, . . .  ,  pm}. 
i:  foreach  vehicle  j  £  {1, . . . ,  m}  \  {i\,  do 
2:  |_  Determine  Uij  using  Algorithm  2. 

3:  Vj  =  Dj=l,.. Uij- 


Proposition  4.2  (Pairwise  dominance  region)  Given 
distinct  locations  p j  =  (A),  Yj),  p j  =  ( Xj,Yj ),  if  a  target 
arrives  at  (x,  0),  where  x  £  Uij  generated  by  Algorithm  2, 
then  T(pi,x)  <  T(pj,x). 

Similar  to  pairwise  dominance  regions,  we  introduce 
the  concept  of  dominance  region  Vi  £  P([0,  W])  for  the 
zth  vehicle,  for  every  i  £  {1, . . .  ,  to},  which  is  the  set  of 
X-coordinates  of  target  locations  for  which  p j  takes  the 
minimum  time  to  intercept  of  all  vehicles. 

Assuming  complete  communication  between  vehicles, 
Algorithm  2  is  extended  to  determine  the  dominance  re¬ 
gion  for  a  vehicle  by  (i)  determining  pairwise  dominance 
regions  between  vehicles  and,  (ii)  taking  intersection  of  all 
pairwise  dominance  regions,  as  presented  in  Algorithm  3. 

The  next  result  for  Algorithm  3  follows  due  to  disjoint 
interiors  of  dominance  regions,  and  due  to  Proposition  4.2. 

Proposition  4.3  (Optimality  of  dominance  regions) 

Given  distinct  vehicle  positions  and  a  target  arrival, 

(i)  the  dominance  regions  generated  by  Algorithm  3  form 
a  partition  of  the  generator. 

(ii)  The  time  taken  to  reach  the  target  is  minimized  by 
the  vehicle  whose  dominance  region  contains  the  tar¬ 
get  arrival  location. 

It  is  possible  for  the  dominance  region  of  a  vehicle  to  be 
empty.  For  instance,  when  one  of  the  vehicles  is  very 
far  from  the  generating  line  (cf.  first  part  of  Figure  8). 
However,  one  condition  under  which  every  vehicle  has  a 
non-empty  dominance  region  is  when  all  vehicles  have  the 
same  F-coordinate.  For  a  general  set  of  locations,  Figure  6 
shows  a  dominance  region  partition  with  three  vehicles. 


Figure  6:  Dominance  region  partition  induced  by  three  vehicles. 


Let  £  :=  [0,  W ]  xR>0,  V([0,  W])  be  the  set  of  all  subsets 
of  [0,  W],  B(r)  be  the  closed  ball  of  radius  r  around  the 
origin  and  +  denote  the  Minkowski  sum  of  two  sets.  The 
domain  of  a  set-valued  map  F  :  X  =1  Z  is  the  set  of  all 
q£l  such  that  F(q)  ^  0.  F  is  said  to  be  upper  (resp. 
lower)  semi-continuous  in  its  domain  if,  for  every  q  in  its 
domain  and  for  every  e  >  0,  there  exists  a  <5  >  0  such 
that  for  every  z  £  q  +  B(6),  F(z)  C  F(q)  +  B(e)  (resp. 

.F(q)  C  F(z)  +  -B(e)).  F  is  continuous  in  its  domain  if  it 
is  both  upper  and  lower  semi-continuous. 

The  pairwise  dominance  region  between  p,  and  p j  is 
a  set  valued  function  Uij  :  £2\Sij  =1  T({0,W])2 ,  where 
Sij  C  £ 2  is  the  set  of  coincident  locations  for  pj  and  p j. 
Similarly,  the  dominance  region  partition  for  vehicle  i  is 
a  set-valued  map  Vj  :  £m  \  Si  =4  V([ 0,  lF])2(m_1),  where 
S  C  £m  is  the  set  of  vehicle  locations  in  which  at  least  one 
other  vehicle  is  coincident  with  pj. 

Proposition  4.4  (Continuity  of  dominance  regions)  (i) 

For  every  distinct  i  and  j  in  the  set  {1, . . . ,  to}  ,  the 
set  valued  map  Uj  is  continuous  in  £2  \  Sij . 

(ii)  For  each  vehicle  i  £  {1, . . . ,  to},  the  set  valued  map 
Vi  is  continuous  on  its  domain. 

Proof:  The  roots  of  Eq.  (7)  which  is  a  quadratic  in 
x,  vary  continuously  with  pj  and  p;.  Thus,  the  map  Uij 
is  continuous  in  £2  \  Sj. 

The  domain  of  Vj  is  contained  in  the  domain  of  Uj  for 
every  j  ^  i.  From  part  (i)  of  this  Proposition,  for  every 
j  ^  i,  the  set-valued  map  Uj  is  upper  semi-continuous  in 
£2 .  Thus,  for  every  j  ^  i,  at  every  q  in  the  domain  of 
Vj  and  for  every  e  >  0,  there  exist  Sj  >  0  such  that  for 
every  z  £  q  +  B(Sij),  Uij( z)  C  Uj( q)  +  B(e).  Given  an 
e  >  0,  by  the  choice  of  <5j  =  min{5jj,  V)  ^  i},  we  obtain 
that  for  every  z  £  q  +  B(8i),  Vj(z)  C  Vj(q)  +  B(e).  Thus 
Vj  is  upper  semi-continuous.  Lower  semi-continuity  of  Vj 
is  established  similarly  and  the  result  follows.  ■ 

4-2.  Minimizing  the  Expected  Constrained  Travel  Time 

For  distinct  vehicle  locations,  Eq.  (3)  can  be  written  as 

m  r. 

Texp(pi,...,pm)=£/  T(pi,x)(j)(x)dx,  (8) 

i= i 

where  Vj  is  the  dominance  region  of  the  zth  vehicle.  The 
gradient  of  Texp  is  computed  using  the  following  formula. 

Lemma  4.5  (Gradient  computation)  For  all  vehicle  con¬ 
figurations  such  that  no  two  vehicles  are  at  coincident  lo¬ 
cations,  the  gradient  of  the  expected  time  with  respect  to 
vehicle  location  p j  is 

dTexp  f  8T 

— - =  /  —{Pi,x)(f{x)dx. 

oVi  Jvt  dpi 

Akin  to  similar  results  in  Cortes  et  al.  (2004,  2005);  Bullo 
et  al.  (2009),  the  proof  involves  writing  the  gradient  as  a 
sum  of  two  contributing  terms.  The  first  is  the  final  ex¬ 
pression,  while  the  second  is  a  number  of  boundary  terms 
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Algorithm  4:  Lloyd  descent  for  vehicle  i 


1: 

2: 

3: 

4: 

5: 

6: 


Assumes:  Distinct  locations  {pi, . . . ,  pm}  G  £ m 
foreach  time  £  G  N  do 

Compute  Vi(£)  by  Algorithm  3  as  a  function  of 

{Pi  (*),••■,  Pm  (*)} 
if  Vi{t)  is  empty,  then 

Move  in  unit  time  to  (A},  4}  —  min{l,  Y)}) 
else 

For  r  £  [t,  f  +  1],  move  according  to 
P i(r)  =  -  sat  f  [  1^-T(pi(T),x)<f>(x)dx' 


which  cancel  out  due  to  continuity  of  T  at  the  boundaries 
of  dominance  regions.  The  complete  proof  is  presented 
in  ?. 

For  z  £  R2,  define  the  function  sat  :  R2  — >  R2  denote 
the  saturation  function,  i.e. ,  if  ||z||  <  1,  then  sat(z)  =  z; 
otherwise,  sat(z)  =  z/||z||.  Inspired  by  the  established 
Lloyd  algorithm  (cf.  Bullo  et  al.  (2009)),  we  present  a 
discrete-time  descent  approach  in  Algorithm  4. 

We  define  the  following  vehicle  configuration. 

Definition  4.6  (Critical  configuration)  A  set  of  loca¬ 
tions  {PD-.Pn,}  is  a  critical  dominance  region  configu¬ 
ration  if  for  all  i  £  {1, . . . ,  to}, 

Pi  =  argmin  /  T(z,x)(f>(x)dx, 
z  e£  J\>i 

where  {Vi, . . . ,Vm}  is  the  dominance  region  partition  in¬ 
duced  by  {pi,...,pm}. 

We  now  state  the  main  result  of  this  section. 

Theorem  4.7  (Convergence  of  Lloyd  descent)  Let  7  : 

N  — >  R2m  be  the  evolution  of  the  m  vehicles  according  to 
Algorithm  f  and  assume  that  no  two  vehicle  locations  be¬ 
come  coincident  in  finite  time  or  asymptotically.  The  fol¬ 
lowing  statements  hold: 

(i)  the  expected  travel  time  t  1— >  Texp(7(t))  is  a 
non-increasing  function  of  time; 

(ii)  if  the  dominance  region  V)  of  any  vehicle  i  is  empty 
at  some  time,  then  V,;  will  be  non-empty  within  a 
finite  time;  and 

(iii)  if  there  exists  a  time  t  such  that  every  dominance  re¬ 
gion  is  non-empty  for  all  times  subsequent  to  t,  then 
the  vehicle  locations  converge  to  the  set  of  critical 
dominance  region  configurations. 

Proof:  We  begin  by  showing  statement  (i).  In  ev¬ 
ery  iteration  of  Algorithm  4,  step  1.  does  not  increase 
the  expected  time  Texp  due  to  the  optimality  of  the  dom¬ 
inance  region  partition,  by  Proposition  4.3.  Step  4:  does 
not  change  the  Texp  as  the  associated  dominance  region  is 
empty.  Finally,  step  6:  does  not  increase  Texp  as  the  vehicle 


is  moving  along  the  gradient  descent  flow  of  Texp.  Thus, 
the  expected  time  is  non-increasing  under  Algorithm  4. 

Statement  (ii)  follows  from  the  fact  that  whenever  Vi  = 
0  for  vehicle  i,  due  to  step  4:,  vehicle  i  reaches  the  generator 
after  finite  time  and  therefore  has  a  non-empty  V*. 

For  non-empty  Vi,  let  A  :  X  x  'P([0,  W\)  — ►  X ,  be  the 
flow  map  of  the  differential  equation  at  step  e.  from  time  t 
to  time  t+1.  For  statement  (iii),  consider  the  discrete-time 
dynamical  system  given  by  the  tuple  (X,Xq,A),  where 
X  =  £m  and  Xq  £  £m  is  the  set  of  initial  vehicle  positions. 

We  now  apply  the  discrete-time  LaSalle  Invariance  Prin¬ 
ciple  (Theorem  1.19  in  Bullo  et  al.  (2009)),  for  which  we 
verify  the  four  assumptions  as  follows. 

1.  Existence  of  a  positively  invariant  set:  At  every 
iteration  of  step  6:,  each  vehicle  follows  saturated  gradi¬ 
ent  descent  of  a  cost  function  belonging  to  the  class  of 
Eq.  (4)  over  its  dominance  region  fixed  for  the  iteration. 
By  Corollary  3.5,  each  vehicle  remains  in  £  throughout  the 
iteration,  and  therefore  at  all  times.  Thus,  the  set  £m  is 
positively  invariant  for  the  system  (X,Xg ,A). 

2.  Existence  of  a  non-increasing  function  along  A:  Texp 
is  non-increasing  along  A ,  by  statement  (i)  of  this  theorem. 

3.  Boundedness  of  all  evolutions  of  (X,  X0,  A)’.  Since 
gradient  descent  keeps  the  X  coordinates  bounded  in  [0,  W], 
it  remains  to  show  that  the  Y -coordinates  of  all  vehicles  re¬ 
main  bounded.  Let  us  suppose  the  contrary.  Then,  there 
are  two  cases:  (a)  at  least  one  vehicle  has  its  location 
bounded  and  at  least  one  other  vehicle,  say  vehicle  k  moves 
so  that  Yk  grows  unbounded;  or  (b)  all  of  the  vehicles  move 
so  that  their  F-coordinates  grow  unbounded.  In  case  (a), 
after  finite  time,  the  dominance  region  Vk  becomes  empty, 
thus  contradicting  the  assumption  of  statement  (iii)  of  this 
theorem.  If  case  (b)  occurs,  then  Texp  grows  unbounded, 
thus  contradicting  statement  (i)  of  this  theorem.  Thus,  all 
evolutions  of  (X,X0 ,A)  are  bounded. 

4.  Continuity  of  Texp  and  A :  Continuity  of  Texp  follows 
from  Eq.s  (2)  and  (8).  To  verify  continuity  of  A,  note 
that  whenever  Vi  is  non-empty,  by  Proposition  4.4,  Vi  is 
continuous  with  respect  to  vehicle  locations.  Thus,  as  long 
as  Vi  is  non-empty,  A  is  continuous  as  the  integrand  is 
continuous  with  respect  to  vehicle  locations,  and  so  is  the 
region  of  integration  V,  . 

By  LaSalle  Invariance  Principle,  the  evolutions  of 
(X,  X0,A)  converge  to  a  set  of  the  form  Te“p(re)n.A4,  where 
k  is  a  real  constant  and  M.  is  the  largest  positively  invari¬ 
ant  set  in  {x  £  X  \  Texp(A(x))  =  Texp(a;)}.  Since  Texp 
remains  constant  under  action  of  A  for  the  set  of  critical 
dominance  region  configurations,  it  is  contained  in  a  set  of 
the  form  Texp(«0  (~l  M.  If  a  set  of  vehicle  positions  is  not 
critical,  then  Texp  strictly  decreases  under  the  action  A , 
and  therefore  the  set  of  vehicle  positions  is  not  contained 
in  a  set  of  T^p(k)  fl  M.  Thus,  the  vehicle  locations  con¬ 
verge  to  the  set  of  critical  dominance  region  configurations. 

■ 

The  next  result  gives  a  simple  condition  to  identify  an 
unstable  critical  dominance  region  configuration,  which  is 
an  unstable  equilibrium  of  Algorithm  4. 
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Figure  7:  Algorithm  4  for  uniform  arrival  density.  The  vehicles  first 
tend  to  an  critical  dominance  region  configuration  (center  figure).  A 
perturbation  to  their  positions  makes  them  move  to  a  stable  config¬ 
uration  (third  figure). 


Figure  8:  Algorithm  4  for  non-uniform  arrival  density  (black  line). 
Initially,  the  blue  vehicle  has  no  dominance  region.  The  vehicles  tend 
to  a  stable  configuration. 


Lemma  4.8  (Disconnected  partitions  are  unstable) 

A  critical  dominance  region  configuration  is  unstable  if 
some  vehicle  has  disconnected  dominance  region  partition. 

The  proof  involves  perturbing  the  position  of  a  vehi¬ 
cle  that  has  a  disconnected  dominance  region  in  a  critical 
dominance  region  configuration,  and  then  showing  that 
the  resulting  gradient  in  the  X  direction  for  that  vehicle 
takes  the  vehicle  away  from  the  equilibrium  configuration. 
The  complete  proof  is  presented  in  ?. 

4-3.  Simulations 

We  now  present  some  simulations  of  Algorithm  4. 
[Examples  of  critical  locations]  We  consider  two  vehicles, 
and  a  uniform  probability  density  of  target  arrival,  i.e. , 
cf>(x)  =  1  jW.  From  initial  locations  such  as  in  the  left¬ 
most  of  Figure  7  wherein  both  vehicles  having  the  same 
A'-coordinate  of  W/2,  but  different  Y-coordinates,  the  ve¬ 
hicles  asymptotically  approach  a  set  of  locations  shown  in 
the  center  figure.  However,  a  small  perturbation  to  the  po¬ 
sitions  leads  the  vehicles  to  positions  in  the  rightmost  fig¬ 
ure.  From  most  initial  conditions,  the  vehicles  converged 
to  a  critical  configuration  as  in  the  rightmost  figure. 
[Non-uniform  probability  distribution]  We  consider  three 
vehicles  and  the  arrival  probability  density  function, 

6(x)  =  iSx/w2’  ifxe[0,w/4], 

nx>  [2/W-8(:r-  W/4)/(3W2),  if  *  €  ]W/4,  W]. 

From  most  initial  conditions,  the  vehicles  converged  to 
a  critical  configuration  as  in  right-most  part  of  Figure  8. 


target  that  arrives  stochastically  on  a  line  segment.  For  a 
single  vehicle,  we  determined  unique  optimal  placements 
when  target  motion  was  either  constrained,  i.e.,  with  fixed 
speed  and  direction,  or  adversarial.  For  the  multiple  vehi¬ 
cle  scenario  and  with  constrained  motion  targets,  we  char¬ 
acterized  conditions  under  which  a  partition  and  gradient 
based  algorithm  that  takes  the  vehicles  asymptotically  to 
a  set  of  critical  locations  of  the  cost  function. 

A  natural  future  direction  is  to  consider  adversarial 
targets  in  the  multiple  vehicles  scenario.  Another  direction 
is  to  consider  stochasticity  in  the  motion  of  the  target. 
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5.  Conclusions  and  Future  Directions 

We  addressed  the  problem  of  optimally  placing  vehi¬ 
cles  having  simple  motion  in  order  to  intercept  a  mobile 
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